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ABSTRACT 

Asteroseismology is a powerful method for determining fundamental properties of stars. We report the first 
application to a metal-poor object, namely the subgiant star v Ind. We measured precise velocities from two sites, 
allowing us to detect oscillations and infer a large frequency separation of /S.v = 24.25 ± 0.25/iHz. Combining 
this value with the location of the star in the H-R diagram and comparing with standard evolutionary models, we 
were able to place constraints on the stellar parameters. In particular, our results indicate that v Ind has a low mass 
(0.85 ± 0.04Mo) and is at least 9 Gyr old. 

Subject headings: stars: individual {v Ind) — stars: oscillations — Sun: helioseismology 



1. INTRODUCTION 

Asteroseismology is a powerful method for determining fun- 
damental properties of stars. This is because oscillation fre- 
quencies give strong constraints on the internal structure that 
are independent of classical observations. Observations of 
solar-like oscillations are accumulating rapidly, and measure- 
ment have recently be reported for several main-sequence and 
subgiant stars, including a Cen A (Bouchy & Carrier 2002; 
Bedding et al. 2004), a Cen B (Carrier & Bourban 2003; Kjeld- 
sen et al. 2005), (3 Vir (Martic et al. 2004; Carrier et al. 2005b), 
/i Ara (Bouchy et al. 2005), HD 49933 (Mosser et al. 2005) and 
■q Boo (Kjeldsen et al. 2003; Carrier et al. 2005a; Guenther et al. 
2005). All of these stars have solar metallicity or greater. As- 
teroseismology is particularly useful for constraining the evolu- 
tionary status of stars with low metallicity (e.g. D' Antona et al. 
2005). Here, we report the first oscillation measurements of a 
metal-poor star. 

We have observed the subgiant v Ind (HR 8515; HD 21 1998; 
HIP 110618), whose iron abundance is only 3% of solar We 
adopt the following stellar parameters: V = 5.28, Tgff = 5300 ± 
lOOK, L = 6.2Lo and [Fe/H] =-1.4±0.1 (Nissen et al. 1997; 
Gratton et al. 2000). Note that the Bright Star Catalogue (Hof- 
fleit 1982) incorrectly lists this star as a binary with spectral 
types A3V:h-F9V. In fact, v Ind is a single star with spectral 
type GO (Lambert & Mc William 1986). 

2. velocity observations and power spectra 

We observed v Ind in August 2002 from two sites. At Sid- 
ing Spring Observatory in Australia we used UCLES (Uni- 
versity College London Echelle Spectrograph) with the 3.9-m 
Anglo-Australian Telescope (AAT). An iodine absorption cell 
was used to provide a stable wavelength reference, with the 
same setup that we have previously used with this spectrograph 
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(Butler et al. 2004). At the European Southern Observatory on 
La Silla in Chile we used the CORALIE spectrograph with the 
1.2-m Swiss telescope. A thorium emission lamp was used to 
provide a stable wavelength reference, and the velocities were 
processed using the method described by Bouchy et al. (2001). 

With UCLES we obtained 680 spectra of v Ind, with typi- 
cal exposure times of 300 s (but sometimes as short as 200 s in 
the best conditions) and a dead time between exposures of 61 s. 
With CORALIE we obtained 521 spectra, with typical exposure 
times of 360 s and a dead time between exposures of 128 s. 

The resulting velocities, with nightly means subtracted, are 
shown in Fig. 1. As can be seen, the weather was good in Aus- 
tralia but poor in Chile (we were allocated seven nights with 
UCLES and 14 with CORALIE). 

Most of the scatter in the velocities, especially for UCLES, 
is due to oscillations. Figure 2 shows a close-up of the first 
night. Oscillations with a period of about 50 min and variable 
amplitude (due to beating between modes) are visible here and 
throughout the time series. We also see good agreement be- 
tween the two instruments, within measurement uncertainties, 
during the overlap (although we should note that the difference 
in absolute velocities is not known and has been adjusted to 
give the best fit). 

Our analysis of these velocities follows the method that we 
developed for a Cen A (Butler et al. 2004) and a Cen B (Kjeld- 
sen et al. 2005). We have used the measurement uncertainties, 
(T,, as weights in calculating the power spectrum (according to 
w, = 1 /erf), but modified some of the weights to account for a 
small fraction of bad data points. In this case, 4 data points from 
UCLES and 24 from CORALIE needed to be down-weighted. 
The power spectra of the individual time series and of the com- 
bination are shown in Fig. 3. The differences between the 
three panels can be attributed to the effects of beating between 
modes. 
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The low metallicity of v Ind means that the Hnes in its spec- 
trum are fewer and weaker than for stars of solar metallicity. 
This, together with its relative faintness {V = 5.3), means that 
the Doppler precision for v Ind is poorer than for other stars 
observed with the same instruments, such as a Cen A and B. 
We measured the average noise in the amplitude spectrum 
of V Ind at frequencies above the stellar signal (0.6-1 mHz) 
to be 17.3cms-i for UCLES, 27.2 cm s"' for CORALIE and 
14.9 cms"' for the combined data. Using these values, we cal- 
culated the noise per minute of observing time to be 5.9 ms"' 
for UCLES and 9.5ms"' for CORALIE, where the difference 
is due to a combination of factors, primarily the telescope aper- 
ture but also including spectrograph design, sky conditions and 
observing duty cycle. 

The inset in each panel of Fig. 3 shows the spectral window 
(the response to a single pure sinusoid). For the single-site data 
we see sidelobes at ±11.6/iHz that are very strong (51% in 
power for UCLES and 57% for CORALIE). These are due to 
daytime gaps in the observing window. When the data are com- 
bined, the sidelobes are drastically reduced (to 16% in power) 
and are also slightly shifted, occurring at ±10.8 /iHz. In the 
cases of a Cen A and B, we generated power spectra in which 
the weights were adjusted on a night-by-night basis in order to 
minimize the sidelobes. We have not done that for v Ind be- 
cause the sidelobes for the two-site data are already quite low. 

3. LARGE FREQUENCY SEPARATION 

Mode frequencies for low-degree p-mode oscillations are ap- 
proximated reasonably well by a regular series of peaks: 

j/„,/ = Aj/(n + i/ + e)-/(;+l)Do- (1) 

Here n (the radial order) and I (the angular degree) are inte- 
gers, /S.V (the large separation) depends on the sound travel time 
across the whole star. Do is sensitive to the sound speed near the 
core and e is sensitive to the surface layers. See Christensen- 
Dalsgaard (2004) for a recent review of the theory of solar-like 
oscillations. 

The large separation, /S.v, is proportional to the square root 
of the mean density of the star. Scaling from the Sun, for which 
/S.V = 135 /iHz, we expect v Ind to have a large separation of 
about 25 /iHz. In order to search for a regular series of peaks 
from which to measure the large separation, we calculated the 
autocorrelation function of the power spectrum. The result is 
shown in Fig. 4 and the peak at 10.8 /iHz (dashed line) is the 
main sidelobe caused by the daily gaps (see above). We inter- 
pret the peak at 24.5 /iHz as due to the large separation, with 
smaller peaks at ±10.8 /iHz from this being due to daily aliases 
(the peak at 13 /iHz also coincides with half the large separa- 
tion). The peak at 5 /iHz is not easily explained by the regular 
p-mode structure and may reflect departures of a few modes 
from the asymptotic relation in equation 1 . 

As a check of this result, we also measured the large separa- 
tion directly from the light curve data using Bayesian methods 
(see Gregory 2005 for a good introduction to Bayesian analy- 
sis). We modelled the light curve as a sum of sinusoids, where 
the number of sinusoids, together with their amplitudes, phases 
and frequencies, were treated as unknowns. An advantage of 
this approach is that the amplitudes and phases can be inte- 
grated out of the problem at the beginning (Bretthorst 1988), 
so that we only need to fit the frequencies and their number To 
determine the large separation, we chose the prior distribution 
for the frequencies to be periodic, with period \/S.v (considered 
unknown). We fixed the central frequency of this periodic comb 



to be the highest peak in the power spectrum (313.14/iHz). 
We then calculated the posterior distribution for Av, which is 
shown in the lower panel of Fig. 4. There is a single strong peak 
with a value of /\v = 24.25 ± 0.25 /iHz, in agreement with the 
peak of the autocorrelation. More details of this method, which 
is similar to the approach used by Bretthorst (2003) and appears 
to be very promising for this type of analysis, will be presented 
separately (B. Brewer et al., in prep.). 

4. CONSTRAINTS ON THE STELLAR PARAMETERS 

What can we learn about v Ind from our measurement of 
/S.vl Figure 5 shows the location of the star in the H-R dia- 
gram, together with some theoretical evolutionary tracks. The 
box, which is the same in each panel, shows the observed po- 
sition of V Ind from classical measurements. The value for T^s 
(5300 ± 100 K) is the mean of published photometric estimates 
(Nissen et al. 1997; Gratton et al. 2000), where we note the 
large uncertainty in the effective temperature scale for metal- 
poor stars. The luminosity is based on the Hipparcos parallax 
(34.6 ±0.6 mas), with bolometric corrections from Alonso et al. 
(1999). Note that the bolometric correction is a function of ef- 
fective temperature, hence the slope of the box. The diagonal 
dashed lines, which are also the same in each panel, are loci 
of constant radius, calculated from L oc R^T^'s' . We can imme- 
diately see that a measurement of the radius from interferom- 
etry would be valuable in constraining the location of the star 
in the H-R diagram, as has already been shown for other stars 
(Kervella et al. 2003; Pijpers et al. 2003; Kervella et al. 2004; 
Thevenin et al. 2005). 

The curved lines in Fig. 5 are evolutionary tracks for a 
range of masses, using model calculations similar to those 
by Christensen-Dalsgaard (1982). We used a metallicity of 
Z = 0.001 (with hydrogen and helium mass fractions of X = 0.75 
and Y = 0.249, respectively) and the three panels differ in the 
adopted value of the mixing-length parameter (a = 1.7, 1.8 
and 1.9), where the solar value is 0:0 = 1.83. The relatively 
rapid evolution in this subgiant phase means each track can be 
described by a single age, as shown in the figure. Finally, the di- 
agonal lines are contours of constant /\v. We calculated these 
from the evolutionary models by scaling from the Sun (since 
/S.V is proportional to the square root of the mean density). 

We see that our measurement of /S.v significantly constrains 
the parameters of v Ind. This is quantified in Fig. 6, in which 
the thin error bars show the range of each parameter based on 
classical measurements alone (L and Teff), while the thick bars 
show the situation after we have added the constraint provided 
by our measurement of /S.v. Including this constraint reduces 
the uncertainty in both effective temperature and radius, for a 
given value of a, by a factor of four 

What can we say about the mass of v Ind? Even with no con- 
straints from seismology, the requirement that v Ind be younger 
than the universe (13.7Gyr; Spergel et al. 2003) sets a lower 
limit of 0.8 1M0. A star of lower mass would not have had 
time to evolve this far. Note that this limit, which we derived 
from the tracks in Fig. 6, is essentially independent of the mix- 
ing length. This is because the value of a has little effect on 
the fusion rate in the core, and hence on the time for a star of 
given mass to leave the main sequence and enter the subgiant 
phase. Meanwhile, an upper limit on the mass is obtained from 
our measurement of /\v, provided we also set a lower limit on 
a (third panel of Fig. 6). Adopting a plausible lower limit of 
a > 1 .7 gives an upper limit for the mass of O.89M0. 
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We can also set interesting limits on the age of v Ind. From 
the bottom panel of Fig. 6, and again setting a > 1 .7 as a plau- 
sible limit, we see that the age must be at least 9 Gyr. This con- 
firms that V Ind is, indeed, very old and must have been formed 
very early in the history of the Galaxy. The final results of our 
analysis are summarized in Table 1, where we list our best esti- 
mates for the parameters of v Ind, assuming that Z = 0.001 and 
a=1.8±0.1. 

With the data currently available, what constraints can we set 
on the mixing length? The dashed line at an age of 13.7 Gyr in 
Fig. 6 indicates the upper limit set by age of the universe. This 
indicates that the mixing length cannot be greater than 2. 1 . We 
expect much stronger constraints on a, and also on the other 
stellar parameters, to come from the individual oscillation fre- 
quencies. The extraction of these frequencies and a comparison 
with theoretical models is deferred to a future paper (F. Carrier 
et al., in prep.). An accurate measurement of the radius using 
interferometry would also be extremely valuable. 

5. OSCILLATION AMPLITUDE 

The amplitudes of individual modes are affected by the 
stochastic nature of the excitation and by the (unknown) value 
of the mode lifetime. To measure the oscillation amplitude of 
V Ind in a way that is independent of these effects, we have fol- 
lowed the method introduced by Kjeldsen et al. (2005). In brief, 
this involves the following steps: (i) smoothing the power spec- 
trum heavily to produce a single hump of excess power that is 
insensitive to the fact that the oscillation spectrum has discrete 
peaks; (ii) converting to power density by multiplying by the ef- 
fective length of the observing run (4.42 d, which we calculated 
from the area under the spectral window in power); (iii) fitting 
and subtracting the background noise; and (vi) multiplying by 
Aj//3.0 and taking the square root, in order to convert to ampli- 
tude per oscillation mode. For more details, see Kjeldsen et al. 
(2005). 

The result is shown in Fig. 7. The peak amplitude per mode 



is 0.95 m s"' , which occurs at a frequency of v^.^^ = 320 /xHz 
(period 52min). This value of ly^ax is consistent with that ex- 
pected from scaling the acoustic cutoff frequency of the Sun 
(Brown et al. 1991; Kjeldsen & Bedding 1995). The observed 
peak amplitude is 4.6 times the solar value, when the latter is 
measured using stellar techniques (Kjeldsen et al. 2005), which 
is substantially less than the value of 7.3 expected from the 
L/M scaling proposed by Kjeldsen & Bedding (1995) but is in 
good agreement with the (L/Mf'-'' scaling suggested for main- 
sequence stars by Samadi et al. (2005). A measurement of the 
mode lifetimes in v Ind would be particularly useful. 

6. CONCLUSIONS 

We have observed solar-like oscillations in the metal-poor 
subgiant star v Ind and measured the large frequency separa- 
tion. We used this, together with the location of the star in the 
H-R diagram and standard evolutionary models, to place con- 
straints on the stellar parameters. Our results, summarized in 
Table 1, confirm that v Ind has a low mass and a large age and 
represent the first application of asteroseismology to a metal- 
poor star. Further constraints on the parameters, particularly 
the mixing length, should come from comparing individual os- 
cillation frequencies with theoretical models. 
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Table 1 
Parameters FOR v Ind (assuming Z = 0.001 and a = 1.8 ±0.1) 



t^v (^Hz) 


24.25 


± 


0.25 


(1.0%) 


Teff (K) 


5291 


± 


34 


(0.64%) 


M{Mq) 


0.847 


± 


0.043 


(5.1%) 


Age (Gyr) 


11.4 


± 


2.4 


(21%) 


L(Lq) 


6.21 


± 


0.23 


(3.7%) 


R{Rq) 


2.97 


± 


0.05 


(1.7%) 


log(g/cm s"^) 


3.421 


± 


0.016 


(3.8% in g) 


angular diameter (mas) 


0.956 


± 


0.023 


(2.4%) 
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Fig. 1. — Time series of velocity measurements otV Ind from tlie UCLES and CORALIE spectrograplis. Tlie mean of each niglit lias been subtracted. 
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methods (see text). 
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Fig. 5. — H-R diagrams in which the box shows the position of u Ind from classical observations and the diagonal dashed lines are loci of constant radius. The 
curved lines are evolutionary tracks for models with Z = 0.001 and a range of masses, with the three panels differing in the value of the mixing-length parameter, ex. 
The diagonal lines are loci of constant /\v, calculated from the mean densities of the models by scaling from the Sun. 
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Fig. 6. — Parameters of u Ind based on observations and models, as a function of tlte mixing-length parameter. The thin error bars show the range of each 
parameter based on classical measurements alone (L and T^n), while the thick bars include the constraint provided by our measurement of Au. The dashed line at 
an age of 13.7 Gyr indicates the upper limit set by age of the universe (Spergel et al. 2003). 
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Fig. 7. — Smoothed oscillation spectrum of u Ind, showing the amplitude per oscillation mode. 
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